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NATIONAL ADVISORY COMMITTEE FOE AERONAUTICS 


TECHNICAL MEMORANDUM No. 1229 


HEAT TRANSMISSION IN THE BOUNDARY LAYER* 
By L. E. Ealikbman 


Up to the present time for the heat transfer along a curved vail in 
a gas flow only such problems have been solved for which the heat trans- 
fer between the wall and the incompressible fluid was considered with 
physical constants that were independent of the temperature (the hydro- 
dynamic theory of heat transfer). On this assumption, valid for gases 
only, for the case of small Mach numbers (the ratio of the velocity of 
the gas to that of sound) and small temperature drops between the flow 
and the wall, the velocity field does not depend on the temperature 
field. 

In 1941 A. A. Dorodnitsyn (reference l) solved the problem on the 
effect of the compressibility of the gas on the boundary layer in the 
absence of heat transfer . In this case the relation between the tempera- 
ture field and the velocity field is given by the conditions of the 
problem (constancy of the total energy). 

In the present paper which deals with the heat transfer between the 
gas and the wall for large temperature drops and large velocities use is 
made of the above-mentioned method of Dorodnitsyn of the introduction of 
a new independent variable, with this difference, however, that the 
relation between the temperature field (that is, density) and the velocity 
field in the general case considered is not assumed given but is deter- 
mined from the solution of the problem. The effect of the compressibility 
arising from the heat transfer is thus taken into account (at the same 
time as the effect of the compressibility at the large velocities). A 
method Is given for determining the coefficients of heat transfer and the 
friction coefficients required in many technical problems for a curved 
wall in a gas flow at large Mach numbers and temperature drops. The 
method proposed is applicable both for Prandtl number P = 1 and 
for P f 1. 

"Gazodimmi oheckaya Teoriya Teploperedachi . " Prikladaaya Matematika 
i Mekhanika, Tern. S, 1946, pp. 449-474. 
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I. FUNDAMENTAL RELATIONS FOR THE LAMINAE AND TUEBULENT BOUNDARY 
LAYER IN A GAS IN TEE PRESENCE OF HEAT TRANSFER 
1. Statement of the Problem 


We consider the flow over an arbitrary contour of the type of a 
wing profile in a steady two-dimensional gas flow (fig. l). For 
supersonic velocities we take into account the existence of an oblique 
density discontinuity (compression shock) starting at the sharp leading 
edge or a curvilinear head wave occurring ahead of the profile. For 
subsonic velocities we assume there are no shock waves (value of the Mach 
number of the approaching flow is less than the critical). 


We denote by u, v the components of the velocity along th9 axes 
x, y, where x is the distance along the aro of the profile from, the 
leading edge, y is the distance along the normal, T is the absolute 
temperature, p the pressure, p the density, p the coefficient of 
viscosity, X the coefficient of heat transfer, e the coefficient of 
turbulence exchange, Xfc the coefficient of turbulent heat conductivity, 
Cp the specific heat, and J the mechanical equivalent of heat. 


T 


T + 


u c 

2Jc, 


is the stagnation temperature 


a 



is the velocity of sound 


is the adiabatic coefficient 


is the Prandtl number. The remaining notation is explained in the text. 
The values of the magnitudes in the undisturbed flew are denoted by the 
subscript oo, the values of the magnitudes at the wall by the subscript 

The problem consists in the solution of the system of equations 
(reference 2) 
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p = pET, n = CT 31 


( 1 . 1 ) 

( 1 . 2 ) 

(1.3) 

(lA) 

(1.5) 


where E is the gas constant, and C and n are constants. 

In the solution of the system equations (1,1) to (1,5) we start 
out from, considerations on the dynamic and the rmal "boundary layer of a 
finite (hut variable ) thickness. The flow outside the dynamic boundary 
layer approaches the ideal (nonviscous) flow, nonvortical in front of 
the shock wave and, in general, vortical behind the wave. Equation (l A) 
shows that the pressure is transmitted across the boundary layer without 
change, that is, the pressure pq(x) and the velocity on the boundary 
of the layer TJ(x) may be considered as given functions of x. Tha 
flow outside the thermal boundary layer we assume to occur without heat 
transfer, that is, outside the thermal layer and on its boun dar y the 
total energy 1 Q is constant: 

u 2 tt2 

Jc p T + TT = Jc p T co + ~ 2 ~ = Jc p T 0 + -7J- = 1 q (1.6) 


Eram this it follows that the stagnation temperature 
thermal boundary layer has a constant value 



T outside the 


( 1 . 7 ) 


Thus the flow outside the thermal boundary layer for small velocities 
is assumed nearly isothermal while for large subsonic velocities isentropic. 
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For supersonic velocities the entropy is constant up to the shock wave 
while after the wave the entropy is constant along each flow line of the 
external flow about the profile but variable from one flow line to the 
next. 


The boundary conditions of the problem are: 

u = v - 0, T* = T w for y=0 (1.8) 

where T w = T^x) is a given function. In the absence of external heat 
transfer (across the wall) we have instead of the second condition of 

equation (1.8) the oondition 

u = U(x) 

T* = T 00 

where 6 y and Zy are the values of y referring, respectively, to the 
boundary of the dynamic and the boundary of the thermal, layer. 

2. Fundamental Expressions for the Temperatures 

For P = 1 equation (1.3) gives the "trivial integral" T* = Constant. 
Taking into account the second condition of equation (1.9), we obtain 

T =* T 00 (l - u 2 ) (u = -~=^ (2.1) 

The temperature of the wall is equal to the temperature of the adiabatic 
stagnation: 


/stA 


\Sy / 


= 0. Further 


y=0 


for y * S y (x) 
for y = Ay(x) 


(1.9) 


= T 00 = 




( 2 . 2 ) 


The integral 
and internal 

was obtained 


(2.1) corresponding to the case of the absence of external 

/prp\ 

heat transfer in the boundary layer, ( — ) =0 for P = 1, 

\£y/y=0 

on the basis of the solution of A. A. Dorodnitsyn (reference l)« 
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For P = 1, U = Constant, and Ty. = Constant, equation (1.3) is 
likewise integrated independent of the solution of the remaining 
equations of tlie system, and gives tlie so-called Stodola— Crocco integral 

T* = au + b (2.3) 

Imposing tlie "boundary conditions we obtain 

T - T oo [1 - a 2 + (.%, - 1 ) ( - $)] (% = (2.k) 

From equation (2.3) ve also obtain 

(t* = T*-GV, t 0 * = T 00 -T w ) (2.5) 

u to 

The integral (2.4) corresponds to the case where there is similarity of 
the velocity field with the field of the stagnation temperature drop 
(equation (2.5)) and was used in the solution of the problem of the 
flow about a flat plate (reference 3) • 


In any more general case U ^ Constant for ( 

ly) ^ 0 or 

\ 

L 

7/ y=0 J 


T^ ^ Constant or P ^ 1 the integral of (l. 3 ) is not known in advance. 

The existence of the trivial integral is not, however, the required 
condition for the solution of the problem and this fact is fundamental 
for what follows. 

Let the function T*(x, y) or u(x, y) be integrals of the system 
(l.l) to (1.5) satisfied by the boundary conditions (equations (l.8) 
and (1.9)). The temperature at an arbitrary point can then be represented 
in the form 

(t* 

T - T oo \t 00 


T = T 00 1 - u‘ 



for 


( 2 . 6 ) 
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As is easily seen, the integrals (2.l) and (2.h) are particular cases of 

t* 

the second form of the relation (2.6) for 5^. = 1 and — — = 

to 

respectively. This relation permits expressing together with the 
temperature also the density and viscosity as a function of the velocity 
and stagnation temperature drop. 


3. Expressions for the Pressure, Density, and Yiscosity 

The pressure p at any point within the boundary layer is determined 
by the equation of Bernoulli 


K 

P = Po = Po^ 1 - u 2 ) K_1 



(3.1) 


where Pq is the pressure on the boundary of the layer (thermal or 
dynamic depending on which of them is thicker), p qq is the pressure 
on adiabatically reducing the velocity to zero in the tube of flow 
passing through the shock wave. The density p at an arbitrary point 
within the boundary layer is determined by the equation of Btate (first 
equation of (1.5)), equations (2.6) and (3.1) 


P = P, 


02 


1-U 2 + (t„. - i)/i * 

V t o 


5 


(1 - D 2 ) 


It 

K-l 


(3.2) 


where p Q2 is the density on adiabatic reduction of the velocity to zero. 
The viscosity jjl is determined by the equation 


n 


P = POO 


1 - u 2 + (T w - l)fl - i--) 


(n S 0.75) 


(3.3) 


where (i 00 corresponds to the temperature Tqo, that is, is obtained on 
the adiabatic reduction of the velocity to zero. 
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For any point ahead, of the shock we have 

K 

p - poiO- - ° s )' t_1 

1 

p = p 01 (l - U 2 )* 1 (3A) 

where and p 0 ^ are, respectively, the pressure and density on 

adiahatically reducing the velocity to zero up to the intersection of the 
streamline with the shock. 

As a sc ale of the velocities it is possible instead of the assumed 
magnitude i/21q to take the critical velocity a* and the local sound 

velocity a. As is known a* = [/(it — l)/(re + 1) / 21 q . Substituting 
■~7 = and — = M we obtain 

a* 1 a 


= / K -l 


U = 


re + 1 


u = 


(re - l)t^/g 
f l + (re - l)#/2 


M = 


\/( is + l)/2 — ( it i— 1 


(2.5) 


k. integral Relations of the Momentum and Energy in New Yarlahles 
From equations (3*l) and (3*2) we obtain 


g = -Pot, _ . -p 0 ro* 


(k.l) 


-y y 

3h order to distinguish various uses of the symbol A, herein, 
subscripts 1 and 2 have been added by the HACA reviewer in the 
translated version. 
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where is the density on the boundary of the layer (thermal or 

dynamics depending on which is the thicker). We represent equations (l.l) 
and (1.2) in the form 


5 (pu2) - h (puv) = p 0 ro* + |:[(p + .) |] 


~ (puU) + (pvU) - puU« = 0 


Subtracting the previous equation from the above we obtain 

^[f § £-§)] + 

+ to'p (^ - ^ + ro*p (1 - j(ii + «) (Jt., 2 ) 


From equation (3*2) we find 


1_ 

K-l 


Pq ** PQ2^ ) 

fo _ 1 = U 2 - u 2 + (^ - 1)(1 - t*/t 0 *) 

P 1 - U 2 


(*.3) 


Integrating equation (^.2) term by term from y = 0 to y = Ay, if A^. > 6 
and to y = 8y if 8 y > Ay (for definiteness we assume that Ay > By.; 
the same result is obtained if we assume 5 y > Ay), making use of the 

relation (h.3), and taking into account the fact that starting from the 
boundary of the dynamic layer the velocity u is constant along x and 
the friction f is equal to zero, we obtain 



NACA TM No. 1229 


9 


* k C p K 1 - §) 37 + wt ( 2 + r^) / 0 Sy p § C 1 - ft ) 47 


+ uu* 




(^A) 


where 


- Rn + «) - M. (g) 

L Sy J y=0 W/w 


(^.5) 


is the frictional stress at the wall. 

We now represent equations (1.2) and. (l«3 J 1 b. the form 


(put 0 *) + (p^o*) - e* 1GT = 0 


s (put * ) + h (p ' rt * ) - pu vr - If + 0 sr] + (I - x ) M? s| 


Subtracting the second equation from the first, integrating the result 
term by term from y = 0 to y = Zy (assuming as above for definiteness 

that Zy > 6y), and remembering that for Zy > 8^. the heat transfer 

q = \ and the friction r = H are equal to zero on the boundary 
of the thermal layer we obtain 




(b.6) 
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where 


= °P 





Is the intensity of the heat transfer at the wall . 


(4.7) 


In solving the problem of the boundary layer without heat transfer 
between the gas and the wall for large velocities and P = 1, 

A. A. Dorodnitsyn introduced the change in variables 


K 


n 




K-l 

(l-U 2 ) dy 


(4.8) 


Noting that the function under the integral sign in equation (4.8) agrees 
with the expression p/pqq for !„. = 1 , we introduce a new independent 
variable of the analogous equation containing in the function under the 
integral sign the expression p/pqq for general case according to 
equation ( 3 . 2 ) 


n = 



1 

1 -U 2 + (^ - 1)(1 - t*/t 0 *) 


K 

(l - U 2 )” -1 dy 


(*.?) 


For T v = 1 


the relation between the coordinates q and y depends on 


the unknown velocity profile. For — rr = ~ equation (4.9) gives the 

+ u 


change in variables applied to the problem of the heat interchange of 
the plate with the gas flow. In this case the relation between q and 
likewise depends only on the velocity profile. 


7 


As is seen from equation (4.9) in the general case the relation between 
the coordinates q and y depends not only on the velocity profile u(x, y), 
but also on the temperature -drop profile t*(x, y) which likewise is 
not initially known but is determined from the solution of the problem. 
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Replac ing tlie density p in equations (4.4) and (4.6) "by its 
expression (3.2) and passing to the variables i, i) ve obtain 




its) jf *('-!) 


dTj 


+ UU* 1 





UCJ'd^-l ) fA/ t*\ 


drj 


— (4.10) 
P 02 


d_ 

dx 






(4.11) 


•where 5 and A are the values of the variable q referring, 
respectively, to the boundary of the dynamic and the thermal layers. 


We denote the thickness of the loss in momentum and the thickness 
of the displacement in the plane xq, respectively, by 



(4.12) 


we introduce the concept of the thickness of the energy loss (in the 
plane xq ) 


9 



dq 


(4.13) 


This magni tude has a clear physical meaning; namely, the magnitude 9 
characterizes the difference between that total energy which the mass of 
the fluid that flows in unit time throug h a given section of the thermal 
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■boundary layer would have if its stagnation temperature were equal to tie 
stagnation temperature of the external flow and the true total energy of 
this mass. The magnitude 9 thus represents in length units referred 
to the temperature to* the "loss" in total 'energy due to the heat 
transfer. For small velocities the concept of the thickness of energy 
loss agrees with the previously introduced concept of the thickness of 
heat-content loss (reference 4). The magnitude 


= 0Hrp 



(4.14) 


may be c all ed the thickness of the thermal mixing. 


With -the aid of the magnitudes defined by equations (4.12), (4.13), 
and (4.14) we represent the obtained integral relations of the momenta 
and energy (equations (4.10) and (4.11)) in the final form 


d0 . U« 
dx u 


E + 2 (H + l) 


Tjg 

1 - U2_ 


d + 


U‘(F W - 1) 
U(1 - U 2 ) 


E t 9 


p u 2 
02 


(4.15) 


d0 U* 
dx U 


9 


V* 

V 


0 = 


%■ 

P02 Uc p t 0* 




(4.16) 


We note that the integration with respect to y (or t]) may be taken 
from 0 to t» so that the relations (4.15) and. (4.l6) are general for 
the theory of the boundary layer of finite thickness and the theory of 
the asymptotic layer. 

For small Mach numbers (the effect of the compressibility due to 
the temperature drop) the relations (4.15) end. (4.l6) assume the form 


dd 21 

dx U 


(H + 2)d + 


HL Af -A 
U V*o 7 


Hrp9 


d0 U* 

32 + IT 



PqU 2 


> 


(4.17) 
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•where t 0 = T Q — T^, T Q and p 0 are, respectively, the temperature and 

density of the isothe rmal flow outside the thermal "boundary layer (it 
follows from, equations (2.6) and (3*2) if we set appr oxima tely IT = 0, 
that is, according to equation: (3.5) M = 0). The pressure distribution 
over the profile is determined by the equation of Bernoulli for an 
incompressible fluid. 1 

As is seen from equations (h.15) and (lf.l6) and also from what follows, 
in the variables x, tj the equations of the system (l.l) to (1.3) a^ 6 
simplified »nrl approach in principle the corresponding equations for the 
Incompressible fluids. For this reason the fundamental methods of the 
theory of the bo undar y layer in an incompressible fluid may be generalized 
to the case of a body in a gas flow with heat interchange. 

We give below the generalization of the method of Pohlhausen for 
the case of the laminar layer and the logarithmic method of Prandtl— 
ir/rm/n for the case of the turbulent layer. The proposed method of the 
solution of the problems connected with heat interchange permits, of 
oourse, generalization of certain other problems in the theory of the 
boundary layer in an incompressible fluid. 

II. LAMINAR BOUNDARY LAYER WITH HEAT INTERCHANGE 
•RHMH ul fliBW TTO GAS AND 'Tint; WALL 2 
5. Transformation of the Differential Equations 


Assuming in equations (l.l) to (1.3) s =» 0 and t = T — T^, 
substituting the values p, p, and |i accord in g to equations (3*1) 
to (3.3), we .transform. these equations to the new independent vari- 
ables x «= "x • and tj, determined according to equation (h.9). ®i© 

equations of transformation of the derivatives will be 

k ' 

s_ = + ^ a_ d (1 - u 2 )*- 1 s_ 

Sx Sx + Bx av ¥ 1 - u 2 + (T* - 1)(1 - t*/b 0 *) 


1 Fram equations (3»l) and (3«5) we have 


J3 2 


poo _p o = sf 1 + 1 (* - 1 ^] K 1 -4 p 0 * : r-( 1 + t^ + £ 2i + •••) 


whence setting M = 0, we obtain 

p 0 + g (pqE 2 ) - Constant 

s In -what follows we restrict ourselves to the case P = 1. 
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We obtain (introducing the notation 



Su 

a 


V 



1 -u 2 + (T^- 1)(1 - t*/t 0 # ) 

UU» 

1 -V s 



If it is assumed approximately that n = 1, the obtained system can be 
still further simplified. 


6. Generalization of the Method of Pohlhausen 

We represent the Telocity profile and the stagnation temperature- 
drop profile by the polynomials 


^^© + ^) 2 + A <") 3 + A < f ) 


K5.1) 


(5.2) 


(5.3) 


( 6 . 1 ) 
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= Ka) + B ^aT +s ll) + Bl (£> 


( 6 . 2 ) 


To determine the coefficients of the polynomials we set up the conditions 


,_J> S + g S. J— *L -1- * + ^ for 1=0 

B( n /s ) 2u A S( n /A) V 5 ^/ 8 ) 11 s 


(6.3) 


where 


_ 2-n 


^2 = 


8*TT% 


v Q2 (i -u 2 )* -1 


fl 


3 a = (l - n) 


^w ~ t qq 
® w 


Condition (6.3) follows from the equation of motion (5.3) since 
u = t = 0 for t] = 0. Further 3 


a - 1, 

IT J 


A_u = o 

U 


a(tj/8) u bctj/s) 


?£=° fOT J- 1 


(6A) 


From equation (5-3) for q = 0, u. = v = 0 we obtain 


> 2 t* 


+ a 


d(q/A) 2 *0* [S(q/A) *0 


B V 


1 2 


J 


= 0 for ^ = 0 


(6.5) 


Similarly to conditions (6.^) for the profile u/U we take for the 
profile t*/to* the conditions 


to* 8(n/A) V 


p ^ 

— “ — onr=° for a = 1 ( 6 *6) 

B(q/A) 2 *0* A 
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By differentiating equations (5»l) and (5.3) with, respect to tj 
with the subsequent equating of tj to zero, it is easy to obtain the 
conditions for the third derivatives of u/fr and t*/tQ* at the wall. 
These conditions may be useful for various aspects of the method of 
Pohlhausen. Using conditions (6.5) and (6.6) we obtain 


B 1 - 


3 -fT= 


12a 


a 


b 2 = 6 - 3B x , b 3 = -6 + 3B 1# B4 = 3 - B-l ( 6.7) 


From conditions (6.3) and (6. if-) we now find 


An = 


.12 + Xr 


1 6 - (3 - i/9 - 12a)8/A 


, A 2 = 6 — 3 Aq, A^ = — 8 + 3 A}_, A4 = 3 ~ A]_ 


( 6 . 8 ) 

In relations (4.15) and (4.l6) there enter, besides d(x) and 0(x ), 
the four unknown functions t , H, q^., %>. In constructing the profiles 

there were also introduced the auxiliary functions 6 and A. The 
required six additional equations are obtained by substituting equa- 
tions (6.1) and (6.2) in equations (4.5), (4.7), and (4.12) to (4.l4). 

We obtain 


6* = -3H = 8 


8 ~il ■d ~ 5Al2 + 12/11 + ^ 


20 


5 


1260 


(6.9) 


K 

,K—1 


\ - 6 (1 - ^ " A 1> 9* - VooA^ 1 <1 - o 2 )* 1 E 1 


8 - Bn 


A* = 6Eqi = A - g0 , 8 = A(M X + JJ 1 A 1 ) 


( 6 . 11 ) 


where for — < 1 
o 


% - todgf - * 3*10, »i - Mtf + + h(d 


6 - B 


”1 “ ~W~’ “ 2 


1 ^ 16 - 3B X 

Do = 


ho = 


5 - ®i „ 24 - 5 B x 


420 * 3 280 


b 4 = 


9 


2520 
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itom equations (6.10 ) and (6.8) it is seen that the point of sepaxation 
of the Tam-tm -r layer for the given "boundary conditions is determined 
by the condition Xg ” 12 * 


7. Determination of the Initial Conditions 

For subsonic flow and also supersonic in those cases where there is 
a head wave in front of the profile a critical point U(o) = 0 is 
formed at the leading edge (x = 0). The latter is a singular point of 
equations (4.15) and (4.l6) in which the derivatives d8/dx, de/dx, 
and so forth, increase to infinity if the initial a, 6 are not 
subjected to certain special conditions. 

Substituting the expressions for t w , q w , A* from equations (6.10 ) 
and (6.11) into equations (4.15) and (4.l6), multiplying the latter 
by 8 a-pd A, respectively, and equating to zero the coefficients of l/U" 
these conditions are obtained in the form 


XgtH + 2) g- + X2 (T w 




8 - B, 


20 



- V-1 = °* 


(7.1) 
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In tlie absence of heat transfer (T w = l) the first relation (7*l) 
is a cubical equation in Xg. Of its three roots (7.052, 17-75, and — JO) 
only the root (^2 ) x _q = 7-052 satisfies the physical conditions. This 

value of Xg i 0 the one generally assumed initial condition in the theory 
of the laminar layer for the equation of Earmimr-Pohlhausen. (The equations 
(7.l) may be conveniently regarded as a system for determining the 
initial values of (^2)^=0 (A/8 ) x= o- We present the results of the 

computation of these values for 0 < (l w .) x _o < 5 (for n = 0.75). 


0 

II 

l&f 

0.05 

0.10 

0.50 

1.00 

1.50 

2.00 

3.00 

4.00 

5.00 

II 

0 

0.3 

0.7 

4. 05 

7.05' 

7.5 

6.0 

3.9 

2.6 

1.9 

^ = 1.12 
C 

1.16 

1.18 

1.23 

1.31 

1.50 

1.87 

2.7*+ 

3.70 

4.86 


The graphs of these results are shown in figure 2. 


8, Method of Successive Approximations 

Simple computation for m ulas can be obtained by generalizing the 
method of H. Lyon (reference 5). At the same time we modify the method 
of Lyon with the object of improving the convergence. 

We multiply (4.15) by 23 and have 


2( . H + 2 + a»r. 


dx 


1 — U 2 


U 1 -t]2 


^ rn H““l 
— 

R 02 U 


(1 i 

1 8 

(8.1) 


where 






E 02 


^00 


and L is a characteristic dimension. Setting k = 2(H + 2) we try 
in the function _k to separate a certain principal part constant for a 
given value of T w ; that is, we set 


k = ci + (k — ci), where Cj = ciCT^) 
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After simple transformations we obtain 


2 s *£ + Hl 

dx u 


°1 + 2(®r - 3 ) vp + + 2(5ff " 2 ]* 2 


= T. n 2 — — (1 - U 2 )* 1 


V 


*®02 




(8.3) 


Assuming over a certain part of the boundary layer xo < x < that 

A* _ 

the ratio — = h is constant with respect to x and equal to the 

mean value for the given segment, we set 3 


= c x + 2(1^ - l)h = c ± + 2(T W - 1) f 


1260(8 - B x ) 


20(-5A 1 2 + 12A ± + l¥+) 


( 8 . 1 +) 


Multiplying equation (8.3) by U° we obtain 


— (^ 2 tJ°) + (c - 2) 
dx 


UU* 
1-U 2 




- n— 2 U 0-1 


T w “ " (1 - U 2 )* - 1 


E, 


02 


1 ’v - x (I) 2 (k - °i> 


Considering this equation as linear in we write its solution in 

the form 4 


3 2 B 0 2 = 


2 L “ z ° 


0+1 T^ 1 * -2 U 3-1 (l - Xj2) 


1 -|+ — 
2 K-l 


4> dx 


3 For ^ Constant there is taken in the exponent c a constant 


mean value of Tv for each section. 


( 8 . 5 ) 
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where 


$ = 





- c-l). 


C = 


73^02^ ( 1 ~ TJ 2 ) 


""I 


( 8 . 6 ) 


X=X r 


where C = 0 for Xq = 0. Substituting under the integral (in the 
function $ and exponent o) certain initial values of (1q)q and (A/ 8 ) q 

we obtain a first approximate function d(x). Ehe arbitrariness in the 
choice of the magnitude c^ may be utilized for improving the convergence. 

For this purpose Cj_ must be chosen such that for each value of T w the 
error due to the assumed (inexact) values of Xq and A /6 is a minimum. 


Since $ itself depends little on A /8 it is sufficient to set up 
the condition of little variation of ® with X,g. Neglecting the relatively 
small dependence of the functions d /8 and k. on X, that is, setting 
in the argument Aj_ on which they depend, Xq = 0, we obtain the 
approximate expression 


^,12 + \ 2 
/w <■ ■ ■ ■ — — — 

6 


2 

(Ax) (|) T - xjft f(k) 
x Xa-oV^o 2 Wx 2 =oL 



In order that <D = Constant the coefficient of Xq must be equal 
to zero, whence we obtain 


c l “ ( k A 

X2=0 6 ( ^/ 5 ) x 2 =o 

Ihls dependence of cp on in a wide interval of change of 

the argument (and practically independent of A/S) is close to a linear 

one. For T w = 0.1 we obtain c-, = 9.35 for ^ = 1 (c-, = 9.12 for 

a A — 6 X 

- = 2 , c^ = 9.7 for — = 0 . 5 ). For T w = 1 we obtain c-j_ = 6.26. 

Bounding off the last value to c 2=6 we apply the linear relation 

°1 = 9.5 - ( 8 . 7 ) 

Finally multiplying the equation of energy (equation (4.17)) in nondimensional 
form term by term by 20 U^ (§ = jj) and integrating as linear we obtain 
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<>*=02 - 


u 2 ^ - 1) 2 


c* + 


r T^ 1 2B 1 (T W - l) 2 U(1 - U 2 ) 

Ji o 




(8.9) 


where 


C* 


'e^R^CTw - I) 2 ] 


±=Xq 


and C* = 0 for £ = 0. 

The computation of the dynamic and thermal layers by equations (8.5) 
and (8.9) can be carried out in this sequence. We consider, together 
with the parameter the analogous composite parameters 


^ 2 Eo2 U,T w 2 n 

2 * = (i-fief 


A^opU'a^ 2-11 
X2 A = (l -U2f ' 


- 0gB O2 TJ, ^ 2 ~ ai 

X2 e" ( 1-02 f 



By equations ( 6 . 9 ) and (6.1l) represented in the form 


A5A-J 2 + 12A 1 + l41+\ 2 

X 2,g - ^2^ 1260 / 

X 2 e = ^2 A ( M 1 + w l A l) 2 



there are constructed for the given value of 1^ auxiliary graphs of 
the dependence of Xq on ^2 8 for various values of A,2 A and f° r the 
dependence of X£ A on Xq q for various values of Xq * The magnitudes 
taken as initial in the computation by formula (8.5) are determined from 



22 


NACA TM No. 1229 


the initial data. If U ^ 0 for xq = 0, there are taken the values 4 (X^q = 0 
and (a/S)q = 1. If tf = 0 far xq = 0, there can he taken as the 
initial values the values of ^ n(i ^2^0 = (A/&)q^(X2)q according 

to figure 2. By the values of the functions d(x), that is, X 2 q (x) of 
the first approximation using the initial value (a/&)q (for the succeeding 
approximations it is convenient to use the initial values of Xg^), ^(x) 

of the first approximation is found with the aid of the graph. Further, 
hy equation (8.9) there is computed 0(x), that is, Xg^Cx) of the first 

approximation, making use of X^x) of the first approx im ation and (a/6) q 
( in the succeeding approximations there is used the function Xg(x) following 
and X 2 a (x) preceding) . From the values of X 2 q (x) and Xg(x) with the aid 

of the graph there is obtained X^Cx) . In those cases where there is a 

considerable change in the ratio A*/i3 (or the function T w (x)) the 
computation must be conducted over segments. The required data for each 
succeeding segment are taken equal to the corresponding values obtained 
at the end of the preceding segment . The local Nusselt number (that is, 
the coefficient of heat transfer q^/tQ* reduced to nondimensional form) 
and the coefficient of friction are found from, the equations 



4 In particular for 1^=1 formula (8.5) with (X^q = 0 gives 

(M.0 < Mo, ) . ' 

\ cr/ 


- 2 1 

5^00 = - -5- 0.47 J X D^d - U 2 ) 2 K_1 dx (c = 6) 

u c (i-ue) 1 '"2 0 

whioh agrees with the equation of L. G. Loltsiansky and A. A. Dorodnitsyn 
for the computation of the .laml mr layer without heat transfer (reference 6). 
In the absence of heat transfer and for smw.l 1 Mach numbers we again 
obtain from this the quadrature 



earlier derived by us on the basis of the method of K^man— Pohlhaus en for 
the laminar layer in an incompressible fluid. (See Tekhnika Yozdushnogo 
Flota, No. 5-6, 1942.) 
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For smal 1 Mach numbers equations ( 8 . 5 ) and ( 8 . 9 ) assume the form 



nr-2 


c-1 


$ dx 




(8.12) 


S®R ■ 


(OV/Tq “ 1) (U/0ce) 2 



Ur-1 


2B n 


V^O 7 u “ A 


dx 


(8.13) 


9. Dependence of the Eeynolds Number Eq 2 on 
the Parameters of the Flow 

Taking account of the fact that according to the equation of state 
p 02 ^02 

we represent the parameter B^o in the form 

poi P 01 ^ 


_ _ ' p 02 
*02 - %L 


where Eq^ = 


/2ioLp 01 


^00 


(9.1) 


The parameter E 0 q is expressed directly in terms of the Eeynolds 


number 


Eeo = 


PcoDpOO 
(J-00 


and the Mach number = — of the approaching 

a co 


flow. From equations (3*3) to (3*5) ve find 


S 01 


- Bo, =- (1 -Uj) 


n— 

T-T 2x K-* 1 


K+l 


XI, 




[1 ♦ I (K - Dm. 2 ] 5 ^ 

r |( K - l)Mco 2 


(9.2) 


1 

2 


Eq 2 = BqI = Rqq. For supersonic velocities the ratio P 02 /P 0 I iB 

found from the condition of a line of a flow passing through an oblique 
shock wave (or a head wave) at the leading edge of the body. Considering 
each surface of the profile separately we denote by (3 q the angle which 
the tangent to the surface of the airfoil at any point makes with the 
direction of the velocity of the undisturbed flow andj by op the angle 
which the normal to the surface of discontinuity makes with the same 
direction. From the equations of the oblique shock wave we obtain 
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•02 


1/2 (it + ljM^cos^ V-l 


P 0 1 \l + 1/2 (k -1)0* oos 2 cp / 


'2k p o 

Moo cos^cp — 

+ 1 


1 

K — l \-^-~ K 

K + l) 


( 9 . 3 ) 


tan (cp + p 0 ) = J- 


1/2(k + l)Moo 2 sin 2cp 


2 1 + 1/2(k - l)Moo 2 cos 2 qj 


( 9 . 4 ) 


In the case of a head -ware In front of the "body the direction of the 
Telocity after the discontinuity (for the flow line at the profile) may 
"be considered to coincide with the direction of the Telocity before the 
discontinuity; in equation (9.3) there is in this case to be substituted 
cp=0. 

10. Boundary Layer in the Flow of a Gas with Axial Symmetry 

For any axial flow about a body of reTolution the integral relations 
of the impulse and energy haTe the following form: 


d f 8 y d T ^ dp a 

— J o pu^ dy-t pur dy = — — 6 y r - ly: 


dx 


( 10 . 1 ) 


~ c p pu(t* - t 0 *)r dy = -q^r 


0 


( 10 . 2 ) 


In these equations the usual simplifications were mads; x is the distance 
along the arc of the meridional section, d the distance along the 
normal to the surface, and r(x) the radius of the cross section of the 
body of rotation (the change of the radius Tector within the boundary 
layei is neglected). The boundary conditions of the problem and also the 
assumptions with regard to the external flow are taken to be the same 
as in section 1. Setting up expressions for T, p, p, and p as in 
sections 2, and 3 and introducing the new independent Tariable T) by 
formula (4.9) we obtain the integral relations of the momenta and energy 
in the variables x, q: 
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&$ XT* 


H + 2 + (H + l) 


U 2 

l-i i 2 


a 


U*(T W - 1) 
TT(l - U2) 



PQO^ 2 


(10.3) 


d 9 U* - r*„ t()l * „ %r 

— + 0 + -=r-0 + v" — 9 = — 7 “* 

as u * V Ooo%V 



— ,and so forth| 


(10.4) 


(L is a characteristic dimension, for example, the length, of the holy 
of rotation. ) Restricting ourselves to the case P = 1 and taking the 
expressions (6.l) and (6.2) we obtain a closed system, of equations 
(10.3), (10.4), (6.9), and (6.11) the solution of which we write (for the 
case where there is no shock wave) in the form 


^\o- 


i-F 


U°(l - U 2 ) 2 r2 *- 


l_2.+_£_ 

+ r d - s 2 ) 2 ic “ 1 35 


*0 


(10.5) 


02RQO = 


IJ2(T w - l) 2 f 2 


C* + f ’Lj ar " 1 2B2P(l - U 2 )^ 1 (T w - l) 2 r 2 | dx 

j_0 ■' 


( 10 . 6 ) 


where 


C = 


1—— 

^OqU^I - U 2 ) 2 r 2 


x=xo 


C* = [^RqoU 2 ^-!)^^ _ 


X=XQ 


The method of computation does not differ from the case of the 
two- dime nsional flow. For the coefficients of the heat transfer and 
friction the equations (8.10) and (8.1l) remain valid. In the case of 
the internal problem (flow in nozzles) the Nusselt number and the 
friction coefficient are determined by the equations 


1*00 = 


^oo^O* 


_ ip n_1 

■ L w 


Bn 


(1 - fi 2 )*- 1 i. 


2r 


w 


%o°=— 55 


2 s; n-1 

w 


Boo 


%,r “(1 - u 2 ) 1 


tt2 \K -1 
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where Xqq is the coefficient of heat conductivity corresponding to 
the temperature Tqq. 

III. TURBULENT BOUNDARY LAYER IN THE PRESENCE OE BEAT 
TRANSFER BETWEEN THE GAS AND WALL 
11. Fundamental Assumptions 


The functions H, Hrp, t , and q v entering equations (4.15) 

and (4.l6) are determined ty equations (4.12), (4.l4), (4.5), and (4.7) 
which express them as functions of S and 9 through the medium of the 
velocity profile and the stagnation temperature-drop profile. The 
present state of the problem of turbulence does not permit representing 
the velocity profile (and also the temperature, profile) by a single 
equation which holds true from the wall to the boundary-layer limit. 

The fundamental dynamic and thermal characteristics of the turbulent 
layer can nevertheless be computed with an aocuracy which is sufficient 
for practical purposes. A fortunate property of equations (4.15) 
and (4.l6) which can be predicted on the basis of the results with 
respect to nonoompressible fluids is that the functions H and Bp 
change very little oyer the length of the turbulent layer and the 
functions t w and q w connected with ■d and 9 by the equations are 
little sensitive to the actual conditiona whioh prevail at a given section 
of the boundary layer. Hence H and Hip (and also magnitudes analogous 
'to them) can be taken as constant over x and the relations between t w 
and -S (the resistance law) and between q^. and 6 (the heat— transfer 
law) can be set up starting from the assumption that the conditions at 
the given seotion of the boundary layer do not differ from the conditions 
on the flat plate. On the basis of the derivation of these supplementary 
equations we assume the simple scheme of Karman according to which the 
section of the boundary layer is divided into a purely turbulent 
"nucleus of the flow" and a "laminar sublayer" in Immediate contact 
with the wall. In the latter the turbulent friction and the temperature 
drop are small by comparison with the molecular. We assume that in the 
turbulent "nucleus" the frictional stress is expressed by the formula 
of Prandtl: 


T 



( 11 . 1 ) 


where l is the length of the mixing path. In other words, in 

equations (l.l) and (1.3) we set e = p Z 2 It follows directly from 

dy 
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2? 


this in view of the fact that the turbulence assumption of Prandtl gives 
= Cp6 that the expression for the heat transfer is 


4 


. dT 
= Xt " 


,2 du dT 

57 5? 


( 11 . 2 ) 


The thickness 5^ of the laminar sublayer of the dynamic "boundary 
layer, equal for P = 1 to the thickness &yi of the thermal., sublayer, 
is determined "by the critical Reynolds number (the Harman criterion) 


= a 2 (aw 11.5) (11.3) 

■where is the velocity on the boundary of the l&minar sublayer, p v 

and are the density and viscosity at the wall. 


12. Derivation of the Resistance Law- 

Assuming that as in the case of the noncompresBible fluid a linear 
variation of the velocity in the laminar sublayer is permissible on 
account of the small thickness, we have 

Ty — Mv _ (12 .1 ) 

5 yZ 

Prom equations (11.3) and (12.1) we obtain 


= a (12.2) 

In equation (12.2) we pass to the variable tj. Near the wall on account 
of the smallness of the terms u^ and t*/tQ* we have 


r d-^) 1 * 

i - u2 + - 1) (I - 

J 0 



K_ 

dT] » T„.(l - V^) 1 ~ K n 


(12.3) 


K 


hence 
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where 6^ * is the thickness of the laminar sublayer for the variable 
Further, 

K 

P w = p 02 ^ Hw = PooV 1 

■‘■w 

Substituting these expressions in equation (12.2) ve obtain 


I* -$J 

T s, ^ 



(12.4) 


For the fundamental parameters of the dynamics layer there is here 
introduced the notation 


** - *«*, £ = ^jg (1 - U e ) 2<K " l) 

V T W/ p 02 J 'W 


(12.5) 


Equation (l2.l) is with the aid of equations (3-2) and (4.9) transformed 
into the form 

3 k 

T ° ‘‘osHt'f ■ ^ a = T ; T j(i-o s ) k - 1 (12. 

[l - l 2 + (ly. — l)(l — t*/t 0 *J 


6) 


Sinoe for sm all tj the terms u 2 and t*/tQ* are smfl.,1.1 and the mixing 
path l = ky (k = 0.391) where the coordinate y ^is expressed according 
to equation (12.3), the "generalized" mixing path l near the wall is 
a linear function of tj: 


K 

y = kn -\ts. (1 - n s ) a( ' K_l) (12.7) 

V /s 

In deriving the resistance law in an incompressible fluid a linear 
mixing-path distribution and a constant frictional stress are assumed 
for the entire section of the boundary layer, from the •w all to the 
outer boundary. Actually the mixing path increases at a considerably 
slower rate than according to the linear law and the friction drops to 
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zero as the outer boundary of the layer Is approached. The assumptions 
made act in opposite directions and lead to a satisfactory relation 
between the parameters and £. 

Carrying over this fundamental Idea of the logarithm! o method into 
our present theory we set t = t in equation (12.7) and assume 

the linear law (equation (12. 7 )) for the entire section of the boundary 
layer. We thus assume that as in the case of the noncompressible fluid 
there will be a mutual compensation of the errors committed in the distri- 
bution of t and 1. Integrating equation (12.6) between T) and 6 
we obtain the approximate velocity profile: 

£ = 1 + — In 2. (12.8) 

u k£ 5 

From equations (l2.l) and (12.4) we obtain the velocity at the boundary 
of the l amin ar sublayer 



a 

z 


The condition of the equality of the velocities of the turbulent and 
laminar flows on the boundary of the sublayer gives 


E 8 - Cil^^ (e 5 = UBEQg, C-l = g e-*** = O.32 (12.10 ) 

Making use of the velocity profile (12.8) in equations (4.12) we obtain 


d_J^/T_2N 6* _ 1 

& " K VI 6 = kg 


( 12 . 11 ) 


Eliminating from equation (12.10) and the first of relations (12.11) the 
auxiliary variable 8, we obtain the resistance law: 


- 0!T v n e t£ (l - ^ 


( 12 . 12 ) 


We obtain incidentally also the approximate expression for the 
parameter H: 

h-2 ! L_ 


■8 1 - 2/kg 


(12.13) 
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13* Derivation of the Heat-Transfer Law 

In this section we shall give a generalization to the case of a 
gas moving with large velocities^ of the heat-transfer law earlier derived 
"by us (reference k) for an incompressible fluid. 

We construct the function 

q* = X = g, + i ru (13.1) 

dy o 

For the turbulent nucleus of the flow we have 


dt dt* ,2 du dt* 

^ * dy * dy dy 

Transforming this equation to the variable q we obtain 

_ *¥2 du dt* *¥2 T? 

p 02 dT l d7 1 1 - h 2 + (5^" - 1)(1 - t*/t 0 * 3 


( 13 . 2 ) 


3 K 

(l - u 2 )* -1 


(13-3) 


Near the wall the function q* behaves like q, ^fchat is, differs little 
from the constant value q w , and the mixing path l depends linearly 

on q according to equation (12. 7). The common mechanism of the transfer 
of heat and the transfer of the momentum in the flows along solid walls 
provides a basis in the derivation of the law of heat transfer for 
assuming as before a constant value q* = q v and the linear law (equa- 
tion (12*7)) for the entire thickness of the thermal layer. Substituting the 
expression for du/dq obtained from equation (12.8) and integrating 
equation (13*3) from q to A we obtain the approximate profile for the 
stagnation temperatures 


t* 

t 0 * 


= 1 


+ 



( 13 ^) 
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From equation (13-1), assuming a linear distribution of idle stagna- 
tion temperatures in the laminar sublayer, we find 

q w = — — (13*5) 

* Ay2 

where t^* is the stagnation temperature at the boundary of the l amin ar 
sublayer . 


As the fundamental therm al characteristics of the boundary layer 
we introduce the following parameters: 


Ervi tt 


Eg = U0 Eq2j 


(1 -u 2 f _1 

T w 00 £ \ 


*00 


ioo t o^ 


(13.6) 


Frcan equations (13 A) and (l3«5) we obtain the stagnation temperature 
on the boundary of the laminar sublayer 



(13.7) 


Equating the stagnation temperatures on the boundary of the iBTninn.r - 
sublayer and making use of condition (12.4) we obtain 


E a = C 1 I w I1 exp(k^)k^ (E a = UABog) (13.8) 

Substituting the expressions for u/U and t*/tQ* according to 
equations (12.8) and (13.4) in equations (4.13) and (4.l4) we obtain 


0 


A 






(13.9) 


From equations (12.10) and (13-8) it follows that ^ = exp(k^ rT , - k£) 
so that we have 5 ± ~ 


0 = _ 1 
A k£ 



(13.10) 
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Eliminating from equations ( 13 . 8 ) and (13.10) the auxiliary par ame ter A 
ve obtain the heat-transfer lav 


E 0 - ClV ^ ^ (3.3.11) 

We obtain incidentally also the approximate expression for the 
parameter H^: 


lM T 

Hrp “ 

iASd - 2M r ) 


( 13 . 12 ) 


14. Solution of the Equation of the Turbulent 
Dynamic Boundary layer 
We represent equation (4.15) in the form 


t U«(H + 1) | - 1) A* 

ax u(i -u 2 ) ^ U(1 -u2) * ^ 


= Il(l 


B 2 )* 1 HR 02 


( 1 ^. 1 ) 


We mate the change in Yariables (reference T){ 


z = e k £(l - 2/k£)k 2 £ 2 (14.2) 

Differentiating this relation with respect to x and using equation (12.12) 
ve obtain 

dx \p« ax ^ ax/ ^ + 2/k 2 £ 2 / 


(ifc.3) 
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From, equations (l4.l) to (A. 3) we obtain 


p K 

dz + gg£ + 1} - + ^ 2 + nE Sd z = 5050 — (1 - U^" 1 (AA) 

^ ud-u 2 ) uCi-u 2 ) 6 ^ ciT w n+1 


The magnitudes H and E which change little along the "boundary 
layer are assumed constant with respect to x. If A /a is considered 
as a known. function of x then equation (A A) is a linear equation with 

respect to z. 

A* 

Assuming a constant mean value of the ratio — = h over a certain 
interval xn < x < xn (in the first approximation we may for the entire 

turbulent layer assume h = = ^5- = E, whioh holds for the plate ), we 

obtain the solution of equation (AA) in the form* 


z = 

where 

C 


/(I -TP? 

rjyiEgC 


01 A /(l- 02 )° 


(A.5) 


* 

= IzT l/(l -U^rJ , c = K(H + 1) + KCT^ - l)h 

L X=Xq 


The constants must be taken equal to 


K = 1.20, E = lA, h = lA, (Ci = 0.326; k = 0.391) 


5 See footnote (3) on page 19 • 



NACA TM Wo. 1229 


3^ 


After computation of the dynamic and thermal layer new mean values 

St E g 

of the variables K, H, h = ■ can he found and, if the deviations 

from the assumed values are considerable , the computation is repeated. 

If, for greater accuracy, the computation is conducted over segments, the 
values of the constants for each succeeding segment are determined by 
the values of R,g, £, Eg, and obtained at the end of the preceding 

segment. In integrating from the point of the transition of the laminar 
into the turbulent state, the magnitude ( z )x=xq is < 3- e ^ erm ^ ne ^ from the 

condition of equality of the initial value of R^ to the value of R^ 
at the end of the laminar segment. In integrating from the leading 
edge, C = 0. By equations (12.2) and (14.2) the auxiliary graphs of the 

functions log (r^^i” 1 ) and log z as functions of can be 

constructed once for all. 

The local friction coefficient is found from the equation 


2t. 




w 


PocPoc 


- S- Fi ♦ I ( « - 1) nil « !!* (* f 4 (1 - 1 ?)*- 1 (i4. 

L e J P 01 VV £2 


6 ) 


For small Mach numbers equation (l4.5) assumes the form 


Z ta 


(Tv/TQ^iu/u, J° 


c + 


KRk 2 

~ 


mrcrf - 


(14.7) 


where 


UooL P n tt 

E = — 0, E 3 - g- *R, C = 


^0 


gffc>‘ 


X=3t- 


15. Solution of the Equation of the Turbulent Thermal Layer 
We represent equation (4.l6) in the form 


age 

dx 


+ 



dtQ* 

dx 


Rg = 





(15.1) 
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Carrying out the change in variables (reference 4) 


Z T = erp 

and making use of the heat— transfer lav ( 13 * 1 ) ve find 


( 15 . 2 ) 


dZm 1 dT n dt n * k2 _ _ x 


K 

K— 1 


( 15 . 3 ) 


= 


i - 1M1 


T 


1 - 2/k^ + 2/k 2 ^ 2 1 


The solution has the form. 


e 


Z T = — 


- if? [ Cx 


c« + 


KrpRogk* 


; f* 5^-i ( 5 v _ x) 1 ® 2_ (1 -B 2 ) 1 ^ 1 as 

'-'Tr\ h b 


(15.^) 


O' . [* T I„ I %V-1 )**]__ 
L J z=lo 


s The equation of energy ( 15 • 1 ) in connection vith the relation 
betveen Eg and ^ in the fom (i 3 .ll) in the case T = Constant is 

an equation vith separable variables so that together vith the solution 
in the form ( 15 A) ve can use its accurate solution 

_ K 

exp(k& r )(kC T - 3) + 2£i(k£r r ) = C + ^ * -2- (1 -U2) K-1 die 

c i^ *o 

vhere 


6 i(h^qi) 



C = 


{exp(k^r p )(k? T - 3) + 2Si(k^rp)l_ 

x=XO 


(for Xq = 0, C = 0) 
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Hie constant must be taken equal to 1.15 (with subsequent check 

by the results of the computation of the thermal layer). ' The magnitude 

___ is determined by equation (l3.ll) from the condition of the 

x=xo 

equality of the value of Eg at the start of the turbulent region to its 
value at the end of the laminar section. In integrating from the leading 
edge, C* = 0. Hie auxiliary graph of the function log zj against k£ip 

can be constructed once for all. 

The local heat transfer is found from the equation 



N = 


*02 U 


^ 


1 + | (k - 1) M* 2 


n 

(1 - u 2 )* -1 


(15.5) 


For small values of" the Mach number equation (15*^) assumes the form 


(T v /T 0 )^(VT 0 - 1)*T 



(15.6) 


1 6 . Determination of the Profile Brag for Subsonic Velocities 

The drag of a wing of infinite span (over unit length of span) is 
obtained from the momentum theorem in the form 


Q - J Pc^JU TO - uj dy = P 00 Uco 2 J ~(l - dT) = P 00 Uoo ^co ( 16 . 1 ) 


where y $ denotes the sum of the momentum-loss thicknesses referred 


to the upper and lower surfaces and computed at a great distance from 
the wing where 5 — and U — ^TIo, For the drag coefficient we have 


2Q 


"xp - 2_ 

pooUoo -L/ 


= Bp. jl + | (K - iJm/J 14-1 


( 16 . 2 ) 
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Ihe problem consists in expressing iS*, in terms of the dyn ami c and 
thermal characteristics of the boundary layer at the trailing edge. 

Since in the wahe behind the body T v = 0 and q^ = 0, equations (4.l6) 

and (i4-.1T ) for the walce assume the form 


dS U* 
dz + TJ 



H + U 2 
i -c 2 


*w-l 

1 - G 2 



0 


( 16 . 3 ) 


d£ 

dx 




0 


We introduce the notation 


(16.4) 


G(i) 


H + U 2 T v - 1 0 „ 


and represent equation (16.3) in the form 


(16.5) 


i ™ . _ (2 + G) i_ J. 

$ dx dx Woo 


(16.6) 


Integrating this equation with respect to x from the trailing edge 
(denoted by the subscript l) to x = 00, we obtain 


■fl TT. rCr 00 tt 

In ~ = (2 + G]_) In =— + / In =7— dG (16.7) 

For Hljss 0, 1^=1 the function G(x) goes over into H(x). For 

an incompressible fluid however the hypothesis of Squire a nd Young 
(reference 8) on the linear character of the dependence ln(TT/CFeo) on H 
holds : 


ln(U /Uj _ ln(Pj/CIco) 
H ^ - E 1 
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Making the analogous assumption 


lnCu/Uaa) ln(Ul/goo) 

G oo “ & G** — 


we obtain 7 


$«, ( G 1 + GA u 

In ■— « 2 + -= ) In 

% \ 2 / tr M 


We write down the expression for Gj_: 


«i - 


^ ^ 2 -??2 A, 


1 - U 2 *1 

It is easily seen that = 1 and as 1, hence 


G„ = 


1 + q» 2 t ^ ~ 1 9 
1 - U 2 + 1 - W 2 * 


(l6.8) 


(16.9J 


(16.10) 


(l6.ll) 


7 Ehe same result can be arrived at from the -following elementary 
considerations. Equation (16.7) may be represented in the form 


In 4 s2 = (2 + G|_) In =!• + (Geo - Gq) In Sa 
*1 ^oa 

where U m is a certain mean value of the velocity U that lies between 
XJj_ and U M . For the usual profile shapes however the ratio Uq/Uco is, 
in general, near unity and since the magnitude 2 + Gq exceeds the 
magnitude Goo — _C>j_ by several times, therefore for any choice of the 
mean value of TJ^ the relative error in the determination of iS^ is 
not large. Taking the geometric mean we again arrive at 

equation (16.9). 
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Replacing in equation (16.9) Gq and Gr M "by their values, we obtain 
finally 



(16.12) 


K, + 5 Hn + 1 

rt = + 


*1* 


Hco 2 


1 - TJ. 


?? 2 


1 -XI 


TT 2 + 2 




1 - U-j 2 


+ t 


1 “IJoo 2 


From equation (10. *»•) it follows that &CJ(T W — l) = Constant, hence 


(v- x) 0 » = (T *i “ 1)e i ^ (l6<13) 

For small Mach numbers equations (l6.2), (l6.12), and (16.13) assume 
the form. 



IT. Boundary Layer in a Gas Flow with Axial Symmetry 

For the turbulent flow about a body of rotation the equations (lO.l) 
and (10.2) in the variables x, y and the transformed equations (10.3) 
and(lO.h) remain valid. Restricting ourselves to the case of the > 

Prandtl number P = 1 and introducing the parameters R^, Rg, and {Yp 

we represent the integral relations in the form. 



bo 
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ggfl + U»(H + 1) + A* _ r» 


^ U(l-U2) 


? (1 - 02) 


dE a _ dt n * n -j 

+ E 0 + _1_ __ E 0 = ™- — HRooCl - U 2 ) 14 - 1 
dx r ix = = 0} -TT 


(17.1) 


(17.2) 


Mating use of the drag lav equation (12.12), the heat— transfer lav 
equation ( 13 .H), effecting the change in variables 


Z = e k f 


{S'* 


J T 




>T 


assuming the little changing magnitudes H, E, constant vith respect 

to x , and also a constant mean value for the ratio — = h, ve obtain 
’ d * 

a system of linear equations the solution of vhich has the form 


/(I - X)2) C 




KRoo* 2 * *****&#& ; * 

c + — „ — r d-ugy 


Jxo i/(l -U 2 ) 0 


tiPvK 1 dx 


(17.3) 


C = 


zU 0 ^ 1 * 3 


L/d - 


Ztj> = 


- 1)^ 


C* + 




'1 J XQ 
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For "fclie frictional stress and the Nusselt number, equations (l4.6) 
and (l5«5) remain in force. In the case of the internal problem 


N 00 “ 




x '00 t 0 


^^ EOOB(1_02/ 


K_ 

K_1 


2t. 


-fOO “ 


v 


Poo 1 


V 2 - 


iLia-ue )*- 1 

£ ®r 


(17.5) 


The theory presented, in particular the integral relations of the 
momenta and energy established in part I, permits determining the 
thermal and dynamic c har acteristics of the boundary layer at a curved 
wall in the most general cases, that is, in the presence of external 
and internal heat interchange. The computation of the boundary layer 
by the equations derived in parts H and 111 on the assumption of the 
Prandtl number P = 1 permits finding directly for arbitrary Mach numbers 
(excluding the interval from to M ro = 1 ) : 

(1) The coefficients of the heat transfer from the wall to the gas 
for a given maintained temperature of the wall through heat supplied 
outside the body and the coefficients of heat transfer from the gas to 
the wall, that is, required for maintaining the heat conduction within 
the body at the given temperature of the wall. 

(2) The distribution of the frictional stress along the wall and 

the profile drag of the wing (in the oase M,,, < for arbitrary 

ratio of the stagnation temperatures and those of the wall. 

For ptmH.n velocities the obtained equations express the dependence 
of the heat transfer and the drag on the ratio of the absolute tempera- 
tures of the flow and the wall (the effect of the compressibility and 
the change of the physical constants due to the heat interchange). 

Tn conclusion we give the results of computation of a single example. 
Tn figure 3 is given the distribution of the velocities of the external 
flow for the super sonio flow about a body with two sharp edges. The 
contour of the body and the position of the discontinuity are also shown. 
The flow was c omp uted by the method of Donov (reference 9)« In figure 4 
are given the curves for the Husselt number H which assure the uniform 
cooling of the surface up to the temperature T^ = 0 . 25 Tqo for 

Ejo = 15 X lO^, = 2, and M* = 6 for the laminar (lower curves) and 
turbulent (upper curves) regimes. 


Translated by S. Reiss 
Rational Advisory Committee 
for Aeronautics 
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